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Abstract 

We investigate nonclassical properties of the field states generated by subtracting any 
number photon from the squeezed thermal state (STS). It is found that the normaliza- 
tion factor of photon-subtracted STS (PSSTS) is a Legendre polynomial of squeezing 
parameter r and average photon number n of thermal state. Expressions of several 
quasi-probability distributions of PSSTS are derived analytically. Furthermore, the 
nonclassicality is discussed in terms of the negativity of Wigner function (WF). It is 
shown that the WF of single PSSTS always has negative values if n < sinh 2 r at the 
phase space center. The decoherence effect on PSSTS is then included by analytically 
deriving the time evolution of WF. The results show that the WF of single PSSTS has 
negative value if 2nt < ln{l — (2n + l)(n — sinh 2 r)/[(29T + l)(ncosh2r + sinh 2 r)]}, 
which is dependent not only on average number *Jt of environment, but also on h and 
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1 Introduction 

Nonclassical Gaussian states play an important role in quantum information processing 
with continuous variables, such as teleportation, dense coding, and quantum cloning. In a 
quantum optics laboratory, Gaussian states have been generated but there is some limita- 
tion in using them for various tasks of quantum information procession pQ. For example, 
when a two-mode squeezed vacuum state (a Gaussian state) with low squeezing is used 
as an entangled resource to realize quantum teleportation, the average fidelity is just more 
(8 ± 2) % than the classical limits. On the other hand, it is possible to generate and manipu- 
late various nonclassical optical fields by subtracting or adding photons from/to traditional 
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quantum states or Gaussian states, which are useful ways to conditionally manipulate non- 
classical state of optical field [2HTT]. Recently, subtracting or adding photon states have 
received more attention from both experimentalists and theoreticians [12H17] . One of rea- 
sons is that photon subtraction can be applied to improve entanglement between Gaussian 
states [I3I19], loophole-free tests of Bell's inequality [20,21], and quantum computing [22] . 
Thus the photon subtraction (a non-Gaussian operation) can satisfy the requirement of 
quantum information protocols for long-distance communication. Nevertheless, the photon 
addition and subtraction have been successfully demonstrated experimentally for probing 
quantum commutation rules by Parigi et al. [23,23]. in fact, they have implemented sim- 
ple alternated sequences of photon creation (addition) and annihilation (subtraction) on a 
thermal field and observed the noncommutativity of the creation and annihilation operators. 

In addition, Olivares et al. [25] theoretically discussed the relation between the photon 
subtracted squeezed vacuum (PSSV), as an output state passing through a beamsplitter, 
and two parameters (the transmissivity of beamsplitter and the photodetection quantum 
efficiency). Then the case of two-mode photon-subtraction is also further discussed in the 
presence of noise [26^127) . Kitagawa et al [28] investigated the degree of entanglement for 
non-Gaussian mixed (pure) states generated by photon subtraction from two-mode squeezed 
vacuum states with on-off photon detectors. For the single PSSV, furthermore, its nonclas- 
sical properties and decoherence was investigated theoretically in two different decoherent 
channels (amplitude decay and phase damping) by Biswas and Agarwal [5]. They indicated 
that the WF losses its non-Gaussian nature and becomes Gaussian at long times in ampli- 
tude decay case. Recently, it is found that consecutive applications of photon subtraction 
(or subtracting a well-defined number of photons) from a squeezed vacuum state result in 
the generation of a squeezed superpositions of coherent state (SSCS) with nearly the per- 
fect fidelity regardless of the number of photons subtracted [B]. The amplitude of the SSCS 
increases as the number of the subtracted photons gets larger. 

It is interesting in noticing that single-mode displaced-squeezed thermal state can be 
considered as a generalized Gaussian state, which has received more attention [29l - f33j . For 
example, phase estimations for squeezed thermal states (STSs) and displaced thermal states 
are presented |29j . which shows that a larger temperature can enhance the estimation fidelity 
for the former. Another example is, for Gaussian squeezed states of light, that a scheme 
is also presented experimentally to measure its squeezing, purity and entanglement [321133] . 
To our knowledge, however, the investigation of photon subtraction from STS (even for 
single photon subtraction case) has not been previously addressed (especially when this 
state interacts with its surrounding environment). In addition, the exact threshold value 
of the decay time has not been explicity given. In this paper, we focus on any number 
photon-subtracted single- mode STS (PSSTS), which is optically produced single-mode non- 
Gaussian states, and explore theoretically its nonclassical properties and decoherence in a 
thermal channel by deriving analytically some expressions, such as normalized constant, 
photon- number distribution and Wigner function (WF). For single PSSTS, it is shown that 
the WF of single PSSTS always has the negative values under the condition of n < sinh 2 r 
at the phase space center (n and r are an average number of thermal state and a squeezing 
parameter, respectively), and that the threshold value of the decay time is dependent not 
only on the average number of environment, but also on n and r. 

In section II, we introduce the single-mode PSSTS, where the normalized factor turns 
out to be a Legendre polynomial with a remarkable result. In Sec. Ill, the nonclassical 
properties of the PSSTS, such as Mandel's Q-parameter, and distribution of photon number 



(related to a Legendre polynomial), are calculated analytically and then be discussed in 
details. In Sec. IV, the explicitly analytical expressions of quasiprobability distributions 
for PSSTS, such as P-distribution, Q-function and WF of the PSSTS, are derived by using 
the Weyl ordered operators' invariance under a similar transformations. Then we derive 
an explicitly analytical expression of time evolution of WF for the arbitrary PSSTS in the 
thermal channel and discuss the loss of nonclassicality in reference of the negativity of WF 
in Sec. V. It is found that the threshold value of decay time corresponding to the transition 
of WF from partial negative to completely positive definite is obtained at the center of the 
phase space, which is not only dependent on the average number 91 of environment, but 
also on n and r. We show that the WF for single PSSTS has always negative value if the 
decay time Kt < \ ln{l — (2n + l)(n — sinh 2 r)/[(29T + l)(n cosh 2r + sinh 2 r)]} (see Eq. (|51|) 
below), where k denotes a dissipative coefficient of interacting with the environment. Sec. 
VI is devoted to calculating the fidelity between the PSSTS and the STS. It is shown that 
the fidelity decreases monotonously with the increment of both photon-subtraction number 
m and the squeezing parameter r. We end with the main conclusions of our work. 



2 Photon-subtraction squeezed thermal state 

At first, let's introduce the photon-subtraction squeezed thermal state (PSSTS). For a 
squeezed thermal field, its density operator is 

Pa = S(r) Pc S\r), (1) 



where S(r) = exp[r(a^ 2 — a 2 )/2] = exp[— \r(QP + PQ)/2] is the squeezing operator [341 
[35] with squeezing parameter r, here the coordinate operators Q = (a + a))/^/2 and the 
momentum operators P = {a — a}) / (\/2i) ( [a, a^] = 1) are introduced as functions of create 
and annihilation operators a) and a, respectively, and p c is a density operator of thermal 
state, 

p c = (l-e°)e^ a , * = ~, (2) 

where k is a Boltzmann constant, and the temperature T is qualified to be a density operator 
of thermal (chaotic) field with trp c = 1. Using the operator identity (3T|,l38j 

= : exp[(e CT - l)cJa\ : = \ exp { J^l (Q 2 + P 2 ) } \ , (3) 

where these two symbols : : and ; ; denote normal ordering and Weyl ordering, respectively, 
and using the Weyl ordering invariance under similarity transformations [37}I38|. which 
means that 

5; (ooo) V 1 = ;5(ooo)5 _1 ; > (4) 

as if the "fence" ; ;did not exist, so S can pass through it, as well as the technique of 
integration within an ordered product of operators (IWOP), one can convert p s to its 
normally ordered Gaussian form [38] (see Appendix A), i.e., 
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t\Ti "~ r \ 2rf 2r 2 



eX P 1-7^2 "7^2 f : > ( 5 ) 



where 



2rf = (2n + l)e 2r + 1, 2r| = (2n + l)e~ Zr + 1, (6) 



n 



which leads to the following relations, 

T 2 - tI = (2n + 1) sinh 2r, (7) 
T \ + r| = (2n + 1) cosh 2r + 1, (8) 
t 2 t 2 2 = n 2 + (2n + 1) cosh 2 r, (9) 

and n = tr (p c a^a) = (e _cr — l) -1 [39] denotes the average photon number of thermal 
(chaotic) field p c in Eq. ([2]). The form in Eq.([5]) is similar to the bivariate normal distri- 
bution in statistics, which is useful for us to further derive the marginal distributions of 

Ps- 

Theoretically, the PSSTS can be obtained by repeatedly operating the photon annihila- 
tion operator a on a squeezed thermal state, so its density operator is given by 

p = C- l a m Ps a^, (10) 

where m is the subtracted photon number (a non- negative integer), and C m is a normalized 
constant with (see Appendix B) 

C m = Tr{a m p s( J m ) = m\D m l 2 P m (b/Vd) , (11) 



which indicates that C m is just related to Legendre polynomial P m (x) (see Appendix B 
(BIO)), and 

B = - [(2n + l)cosh2r - 1] , (12) 
D = n 2 — (2n + 1) sinh 2 r. (13) 

It is noted that, for the case of no-photon-subtraction with m = 0, Cq = 1 as expected. 
Under the case of m-photon-subtraction thermal state with B = n, D = n 2 , and P m (1) = 1, 
C m = m\n m . The same result as Eq.(24) can be found in Ref. |40| . 

Here we should point out that, as Agarwal et al introduced the excitations on a coherent 
state by repeated application of the photon creation operator on the coherent state |41| . 
we introduce theoretically the PSSTS (|10p . In realistic situations, one the other hand, the 
photon subtraction would be done by on/off detector and the tapping beam splitters with 
a non-unity transmittance, which leads to a generated mixed state. For various schemes for 
generating photon subtraction, one can refer to Refs. [Tl l28ll4*2] . 

3 Nonclassical properties of PSSTS 

3.1 Mandel's Q-parameter 

The analytical expression of C m is of importance for further investigating the properties of 
PSSTS. For instance, one can easily calculate 

ata) = Tr(C- 1 a-+ 1 p s a tm+1 ) = % ±i , (14) 
atV) = lr{C-^ 2 p s a^ 2 ) = %±?-, (15) 
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Figure 1: (Color online) The Q-parameter as the function of squeezing parameter r for different 
m = 0,1,2,3,4, 19,20. 

thus the Mandel's Q-parameter is given by 



which measures the deviation of the variance of the photon number distribution of the field 
state under consideration from the Poissonian distribution of the coherent state. If Qm = 
we say the field has Poissonian photon statistics, while for Qm > (Qm < 0)> the field has 
super-(sub-) Poissonian photon statistics. It is well-known that the negativity of the Qm- 
parameter refers to sub-Poissonian statistics of the state. But a state may be nonclassical 
even though Qm is positive as pointed out in Ref. [35] . This case is true for the present 
state. In fact, if Qm is positive, it does not mean that the state is classical. In such cases, 
we have to use other parameters to test the non-classicality |43| . From Fig.l, one can see 
clearly that for odd number m, Qm becomes negative when the squeezing parameter r is less 
than a certain threshold value which decreases as m increases. Differently from the case of 
odd number m, Qm is always positive for even number m. It is necessary to emphasize that 
the Wigner function (WF) has negative region for all r, and thus the PSSTS is nonclassical. 
In addition, when the average photon number n is larger than a certain threshold value, 
Qm is also always positive. Without loss of generality, thus, we consider only the (ideal) 
PSSTS in a thermal channel in our following work. 

3.2 Photon-number distribution of PSSTS 

Next we discuss the photon-number distribution (PND) of PSSTS. Noticing a) m \n) = 
W (m + n)\/n\ \m + n) and using the un-normalized coherent state |a) = exp[aa^] |0), |44|I45| 
leading to |n) = -rr^r \ot) | a =o , ({ft |«) = e a/3 ), as well as the normal ordering form of p s 
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in Eq. (|1J) 

V(n) 



the probability of finding n photons in the field is given by 

= (n\ p\n) = C" 1 (n| a m p s a) m \n) 
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where Ai and A 2 are defined by 
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In a similar way to deriving Eq.(Bll), we finally obtain 
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where -P m + n (2?) is Legendre polynomial in (BIO), and 
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In particular, when m = (Co = 1), Eq. (|20j) reduces to 
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which is just the PND of STS which seems a new result; while for r = (t\t 2 = n + 1, A\ 



n/(n + l), E = n 2 /(h + 1) 2 , P m (l) 
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Eq. ([20j) becomes 



7) 



m+n 



which is the PND of m— photon-subtracted thermal state which also seems a new result, 
and the PND (n n j (n + l) n+1 ) of thermal state without photon-subtraction [46|, l47j . 

In Fig. 2, the PND is shown for different values (n; r) and m, from which we can see 
that by subtracting photons, we have been able to move the peak from zero photons to 
nonzero photons (see Fig. 2 (a)-(c)). The position of peak depends on how many photons 
are annihilated and how much the state is squeezed initially. In addition, the PND mainly 
shifts to the bigger number states and becomes more "flat" and "wide" with the increasing 
parameter r and the average photon-number of thermal field p c (see Fig. 2 (b) and (d)). 




4 Quasiprobability distributions of PSSTS 

In this section, several quasiprobability distributions of PSSTS are derived in order to 
provide a convenient way for studying the nonclassical properties of fields. 



4.1 P distribution 

We first calculate the Glauber-Sudarshan P distribution function [38] of PSSTS. For this 
purpose, we start from the anti- normal ordering form of p in Eq. (llOp . Recalling the 
integration formula converting an operator O into its anti-normal ordering form [49j with 

anti- normal ordering : :, i.e., 

= J^-': <-/3| 6 |/3) exp (|/?| 2 + p*a - /3a f + a f a) :, (24) 

where \(3) is a coherent state, one can obtain the anti-normal ordering form of the squeezed 
thermal state p s by substituting Eq.([5|) into Eq. (|24p and using the integration formula (B7). 
Such as 
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thus the anti-normal ordering form of p in Eq. (]10p is 



p= -^;a m exp 
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which leads to the P-function P (a) of PSSTS, 

P(a) = C- 1 \a\ 2m P (a), 
where Pq ( a ) is the P-function of STS, 

1 



P (a) 
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(26) 



(27) 



(28) 



It is interesting in noticing that when r = 0, Eq. (l28h becomes P (a) = \a\ 2m e~' Q ' ^ n /(nC m ), 
which is just the P-function of m— photon-subtraction thermal state which seems a new 
result. From Eq. (|27p one can see that the P-representation of density operator p can be 
expanded as 



C£j *JL\z\*»Pb{z)\z)(z\ 



(29) 



which is a non-Gaussian function due to the presence of \z\ 2m . 



4.2 Q-function 

The Q-function is the absolute magnitude squared of the projection of a state of the field 
onto a coherent state (a\, defined by 



Q (a, a*) = — (a I p la) . 

Substituting Eq. ([29|) into (p0|) . we can obtain 

Q (a, a*) = R m (a, a*) Qo (a, a*) , 
where Qq (a, a*) is the Q-function of STS, 
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which seems a new result not reported before, and R m (a, a*) is a factor generated from 
the photon-subtraction, i.e., 



R m (a, a*) = C fn 1 y] 
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where M = [(2n+l) sinh 2r]/[4 (n 2 + (2n + 1) cosh 2 r)], and O = 2n(n+l)/[(2n + 1) sinh2r]. 
Eq. (l3ip indicates that the Q-function of PSSTS is also a non-Gaussian type due to the pres- 
ence of R m (a, a*) and always positive since O > 0. In particular, when m = 0, R m (a, a*) = 
1, thus Q(a,a*) = Qq (a, a*), as expected. 



4.3 Wigner function 



Next, the P-function is applied to deduce the WF of PSSTS. The partial negativity of WF 
is indeed a good indication of the highly nonclassical character of the state. Therefore it is 
worth of obtaining the WF for any states. The WF W (a, a*) associated with a quantum 
state can be derived as follows [361 : 



(34) 



W (a, a*) = tr[pA (a, a*)], a = (q + ip) /V2, 
where A (a, a*) is Wigner operator, whose coherent state representation is 



A (a, a* 
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(35) 



where \/3) = exp(— |/3| 2 /2 + fiefi) |0) is the coherent state. Using the vacuum projector 
|0) (0| = : e~ a a : , and the IWOP technique [37] one can put Eq. (135f) into its normal ordering 
form, 

' f -2 (a) - a*) (a - a) 



A (a, a* 



— : exp 
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Thus substituting Eqs.flHD, ® and O into Eq.flMD, we can finally obtain the WF of 
PSSTS (see Appendix C), 
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where W (a, a*) is the WF of STS, 
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where /3 = [2a*(n — sinh 2 r) + asinh 2r]/{z[(2n + 1) sinh 2r] 1//2 }. Eq. (|37p is the analytical ex- 
pression of WF for PSSTS, related to single- variable Hermite polynomials. It is obvious that 
there does not exist negative region for WF in phase space when n > sinh 2 r which is agree- 
ment with Eq.(28) in Ref. [50]. In particular, when m = 0, Fq (a, a*) = 1, Eq. (I37p becomes 



W (a, a*) = W (a,a*); while for r = 0, note C n 
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which corresponds to the WF of m-photon subtracted thermal state [40] , and can be checked 
directly from Eq.(C6). In addition, for m = 1, (single-photon-subtracted squeezed thermal 
state (SPSSTS)), C x = B $2$, the special WF of SPSSTS is 



Wt [a, a*) = Fi (a, a*) Wo (a, a*) , 



(41) 





Figure 3: (Color online) Wigner function distributions W (a, a*) of PASSTS for different (n, r) and 
m values (a) n = 0.1, r = 0.5, m = l;(b) n = 0.1, r = 0.5, m = 2;(c) n = 0.2, r = 0.5, m = l;(d) 
n = 0.1, r = 0.8, m = 1. 



where 

Fl ( °' °* } = (2n+l)f? ( (2 ™ + X ) I"!' + ^ " sinh2 r ) 

1 2a* f n — sinh 2 r) + a sinh 2r 1 2 n — sinh 2 r , > 
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Noting B > 0, thus from Eq. (|4ip one can see that when the factor F\ (a, a*) < 0, the 
WF of SPSSTS has its negative distribution in phase space. This indicates that the WF of 
SPSSTS always has the negative values under the condition: n < sinh 2 r at the phase space 
center a = 0, which is similar to the case of single-photon-subtracted squeezed vacuum 07]. 

Using Eq. ([37|l . the WFs of PSSTS are depicted in Fig. 3 for several different values of n, r 
and m in phase space. It is easy to see that the the WF is non-Gaussian in phase space. 
As an evidence of nonclassicality of the state, squeezing in one of the quadratures is clear 
in the plots (see Figs. 3(a) and 3(d)). In addition, we can clearly see that there is some 
negative region of WF, which is another evidence of nonclassicality of the state, and that 
the negative region of WF gradually disappears as the n (or the temperature) increases for 
given r and m (see Fig. 3(a) and (c)). Furthermore, for a larger squeezing, the WF shows 
a smaller minimum negative value at the center of phase space (Figs. 3(a) and 3(d)). For 
two-photon subtracted case, the WF presents two positive peaks and two negative peaks, 
different from the case of single-photon subtracted case. 



5 Decoherence of PSSTS in thermal environment 

When the m-PSSTS evolves in the thermal channel, the evolution of the density matrix can 
be described by master equation [51] 

= k (9t + 1) i/lapa} — a) ap — pa) a) + nVl (2a) pa — aa) p — paa)\ , (43) 



1 n 



where k represents the dissipative coefficient and 9T denotes the average thermal photon 
number of the environment. When 91 = 0, Eq. (|43p reduces to the master equation describing 
the photon-loss channel. The evolution of the WF is governed by the following integration 
equation [52J, 
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where W(a,a*,0) is the WF of the initial state, and T = 1 — e~ . Eq. (jl4"|) is just 
the evolution formula of WF in thermal channel. Thus the WF at evolving time may be 
obtained by performing the integration with an initial value. 

Substituting Eqs. ([37]) - ([39j) into ([4"1"]) . and using Eq.(B7) we finally obtain the evolution 
of WF for PSSTS in thermal environment (see Appendix D, in a similar way to deriving 

Eq.flll])), 

w((,c,t) = F rn ((,c,t)Wo(c,e,t), 



where Wq (£, £*, t) is the WF of squeezed thermal state in thermal channel, 
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(49) 



It is noted that, at the initial time (t = 0), (20t + 1) TVG -)• 1,A 2 -)• 2g ,^ -)• 
(/ 2 , Ai -)■ |52,X ->• 9l) ^ -> 2 5'iC + S^C*, Eqs.flM]) and (g7]) just reduce to Eqs.([3H]) and 
(|39|) . respectively, i.e., the WF of PSSTS. In addition, for the case of m = 1, corresponding 
to the case of SPSSTS, Eq. fiTl) just becomes 



F 1 (c,C,t) = c^ L M' + x 



(50) 



from which one can see that when the factor F\ (£, £*, t) < 0, the WF of SPSSTS in thermal 
channel has its negative distribution in phase space. At the phase space center ( = 0, 
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Figure 4: (Color online) Wigner function W (a,a*,t)oi SPSSTS for r = 0.3, n = 0.05 and 
different rand Kt: (a) Kt = 0.05;(b) Kt = 0.1;(c) Kt = 0.2; (d) ret = 0.5. 



the WF of SPSSTS always has the negative values when x < 0, leading to the following 
condition: 

1 

1 



Kt < Kt e = — In 
2 



2n + 1 n — sinh 2 r 



(51) 

291 + 1 n cosh 2r + sinh 2 r J ' 
which implies that the threshold value nt c is dependent not only on the average number 
91 of environment, but also on the average number n of thermal state and the squeezing 
parameter r (a result different from other discussions about the threshold value nt c in 
thermal channel 0(7] )• The WF of SPSSTS is always positive in the whole phase space 
when Kt exceeds the threshold value Kt c . Actually, Eq. (]5ip is also true for the case with 
any number (m) photon-subtraction (see Eq. (|47p ). From Eq . (|51|) one can clarify how the 
thermal noise (n, 91) shortens the threshold value of the decay time. 

Using Eq. (|45|) we present the time-evolution of WF at different times scales in Fig. 4. 
From Fig. 4, one can see clearly that the partial negative region of WF gradually diminishes. 
At long times Kt — > oo, one has oj — > 0, x ~ > n cosh 2r + sinh 2 r, Ai — > \ (2n + 1) sinh 2r and 
H m (0) = (-l) fe m\5 m: 2k/kl. Thus 

w ^'^ = VWTTf^' (52) 

which is independent of photon-subtraction number m and corresponds to thermal states 
with mean thermal photon number 91. This implies that the system reduces to thermal state 
after a long time interaction with the environment. Eq. (|52p denotes a Gaussian distribution. 
Thus the thermal noise causes the absence of the partial negative of the WF if the decay 
time Kt exceeds a threshold value. In addition, from Fig. 4, it is found that the SPSSTS is 
similar to a Schrodinger cat state. 

6 Fidelity as a non- Gauss ianity measure for PSSTS 

Recently, some quantitative measures to assess non-Gaussianity are proposed |53|,I54]. A 
non-Gaussianity measure may serve as a guideline to quantify the non-Gaussian states. 



Therefore, it is of interest to evaluate the degree of the resulting non-Gaussianity and 
assess this operation as a resource to obtain non-Gaussian states starting from Gaussian 
ones. Here, we examine the fidelity between the PSSTS p s and the STS p. Since the 
STS can be considered as a generalized Gaussian state, the fidelity may be seen as a non- 
Gaussianity measure able to quantify the non-Gaussian character of a quantum state. In 
order to quantify the non-Gaussian character of the PSSTS, we introduce the fidelity by 
defining 

3 = tr(p s p)/tr(p 2 ), (53) 

where p s and p are the squeezed thermal state (a generalized Gaussian state) and the PSSTS, 
respectively. Obviously, when photon-subtraction number m = 0, leading to p = p s , then 
5 = 1 which means that p is a Gaussian state described by p s . 
Using Eqs.([T]) and ([2]), one has 

tr(^)=tr(, c 2 ) = ^L_. (54) 

On the other hand, the fidelity (tr (p s p)) can then be calculated as the overlap between the 
two WFs: 

tr(p s p)=47T J d 2 aWo(a,a*)W p (a,a*) , (55) 

where Wo (a, a*) is the WF of squeezed thermal state p s . Using Eq. ([37|) we may express 
Eq.([55]) as 

tr(p s p) = 4vr J F m (a,a*)W^ (a,a*)d 2 a. (56) 

Then employing Eqs, (j38j) and (C6), similarly to Eq. ffTTj) . Eq. (|56j) may rewritten as (see 
Appendix E) 



m/2 



m\B. 



tr ^' = (2a + l) C „ P "' B '^ B2 '' < 57 > 



where P m (x) is the Legendre polynomial with 

^2/ = , r N2 



„ n (n + 1) „ n (n + 11 9 , 9 , . 

Bi = — i '- cosh 2r, B 2 = — £- - sinh 2 r cosh 2 r. (58) 

2n+l (2n + l) 2 V ' 



Thus the fidelity (J53D for the PSSTS is given by 



C m V ' \DJ p m [B/VD) 



which is an analytical expression for the fidelity between PSSTS and SSTS. We see that 
when m = (the case of no photon-subtraction), $ = 1; while for m = 1 (the case of 
SPSSTS), Eq.(j59|) reduces to 

_ n (n + 1) cosh 2r 

(60) 



(2re + 1) (sinh r + n cosh 2r 



In Fig.5, we plot the fidelity between PSSTS and STS as the function of squeezing 
parameter r for different photon-subtraction number m. From Fig. 5 one can see that the 
fidelity decreases monotonously with the increment of both photon-subtraction number m 
and the squeezing parameter r, as expected. 
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Figure 5: (Color online) The fidelity between PSSTS and squeezed thermal state as the function of 
squeezing parameter r for different photon-subtraction number m = 0,1, 2, 3, 4, 19, 20. (n = 0.2). 
The cases of m =19 and 20 are not identical, but they are almost overlap each other, which can not 
be seen clearly from figure due to the use of thick style for line. 

7 Conclusions and Remarks 

In summary, we investigate the nonclassicality photon-subtracted squeezed thermal state 
(PSSTS) and its decoherence in thermal channel with average thermal photon number 91 
and dissipative coefficient k. For arbitrary number PSSTS, we have, for the first time, 
obtained an analytical express for the normalization factor, which turns out to be a Leg- 
endre polynomial of squeezing parameter r and average photon number n of thermal state, 
a remarkable result. Based on Legendre polynomials' behavior the nonclassical properties 
of the field, such as Mandel's Q-parameter and photon number distribution, are also de- 
rived analytically, Furthermore, the nonclassicality of PSSTS is discussed in terms of the 
negativity of WF after deriving the explicit expression of WF, which implies the highly 
nonclassical properties of quantum states. It is shown that the WF of single PSSTS al- 
ways has negative values if n < sinh 2 r at the phase space center. Then the decoherence 
of PSSTS in thermal channel is also demonstrated according to the compact expression 
for the WF. It is found that the threshold value of the decay time corresponding to the 
transition of the WF from partial negative to completely positive definite is obtained at 
the center of the phase space, which is dependent not only on the average number 91 of 
environment, but also on the average number n of thermal state and the squeezing param- 
eter r. We show that the WF for single PSSTS has always negative value if the decay time 
Kb < 5 ln{l - (2n + l)(n - sinh 2 r)/[(29t + l)(ncosh2r + sinh 2 r)]}. A non-Gaussianity 
measure may serve as a guideline to quantify them for the class of non-Gaussian states, 
where the fidelity decreases monotonously with the increment of both photon-subtraction 
number m and the squeezing parameter r. 

In addition, Mandel's Q parameter does not always indicate a negative value for non- 
classical state. In fact, for the photon-subtracted squeezed states by even number, this 
parameter is positive. Thus the negativity of Q parameter is a sufficient condition to 
distinguish non-classical state from classical one. While for photon subtracted squeezed 
states by odd number, the negativity of single photon subtracted case is noticeable. To 



1 A 



compare further non-classicality of quantum states for different number subtracted case, 
the measures based on the volume of the negative part of the Wigner function [55], on the 
nonclassical depth |56| and on the entanglement potential |57] may be other alternative 
methods. Non-classical state introduced in this work will maybe used in combination with 
other non-classical states such as entangled states. 

On the other hand, we should mention that for a photon-subtracted squeezed state gen- 
erated with some realistic probability, its non-classicality, in particular, its non-Gaussianity 
would not be always superior to the input Gaussian state. For example, photon-subtracted 
two-mode squeezed vacuum state has more entanglement than initial two-mode squeezed 
state in not so strong squeezing parameter; while for strong squeezing region, its superior- 
ity disappears [28]. Entanglement evaluation investigation for photon-subtracted two-mode 
squeezed thermal state is a future problem. 
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APPENDIX A: Derivation of Eq.© 

Using the operator identity ([3]) and noticing the single-mode squeezing operator yields 
the transformations, 



S (r) Q5 f (r) = e~ r Q, S (r) P5 f (r) = e r P, 



(Al) 



one has 

Ps = (1 - e°)S{r)e° a ' a S\r) 



2(1 



e a + 



(A2) 



which is still in Weyl ordering, in deriving (A2) we have used the Weyl ordering invariance 
under similarity transformations (j3|). According to the definition of Weyl correspondence 
rule [58], i.e., the classical Weyl function f (q,p) of operator p s can be given by replacing 
the Q and by q and p in its Weyl ordered form, respectively, 



ft ^ 2(1 -e"j \< 
J(q,P) = — | -, exp 



e a + 1 



1 



e a + 1 



' -2r 2 , 2r 2~\ 

e q + e p 



(A3) 



then using the relation between p s and Wigner operator A(q,p) , i.e., operator p s can be 
expanded in terms of A (q,p), 



dqdpf (q,p) A(q,p) , 



(A4) 



where the normal ordering form of A (q,p) is given by 

1 



one can see that 



A (<?,£>) = - : exp 

7T 



{q-Qf-{p-Pf 



2(1 - e a ) 
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dqdp exp 
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x : exp 



(q ~ Q) -ip-P) 



e a + 1 
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' —2r 2 , 2r 2 

e q + e p 



>} 



= Eq.©. 



(A6) 



thus we complete the proof of Eq. ([5]) . 

APPENDIX B: Deduction of Eq.flTTD 

Using the completeness relation and p' s s normal ordering form in ([5]), as well as the 
overlap of coherent state, 
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(Bl) 



we have 
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(B2) 

(B3) 
(B4) 
(B5) 

(B6) 



and using the integration formula |59| 
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whose convergent condition is Re(£ ± / ± g) < and Re ^ ^±f±g j < an d noting that 
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one rewritten Eq. (B2) as 

c m = (-Ay 



& 



2 m 



dk m ds r 



■ exp 



.jfe2 _ s 2 _ 

A 



s=k=0 



(_A) m 2 m m! 



m/2] , /_R\m— 2n 

ml 1 1 



2A . 



~ 2 2n (nly (m - 2n)\ 



(B9) 



Recalling the newly found expression of Legendre polynomial (its equivalence to the well- 
known Legendre polynomial's (P m (x)) expression is [60J 
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we derive the compact form for C m , 
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(Bll) 
where 

D = B 2 — 4A 2 = n 2 — (2n + 1) sinh 2 r. (B12) 

Eq.(Bll) indicates that the normalization factor C m is just related to Legendre polynomial. 

Combining Eqs.(B8) and (BIO), On the other hand, one can derive a new formula for 
Legendre polynomial, i.e., 

d 2m ( , 2xrt \ 2 m m\ 



dt m dT r 



exp -r - -r + 



t lT =o (x 2 - 1) 



m/2 



P m (x) . (B13) 



APPENDIX C: Derivation of WF flTTD for PSSTS 

According to Eqs.flMJ), (JHHJ) and ([29]), we have 

W(a,a*) = C m Hr[J ^-\z\ 2m P {z)\z) (z\A(a,a*)} 
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which leads to 
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Then using the integration formula (B7), we can write Eq.(Cl) as following form, 

W ^ = ^^takZt- 6XP ^ (2a + k) ^ + *) 

4 V /Jfc=t=o 
= F m (a, a*) Wo (a, a*) , 

where Wq (a, a*) is the WF of squeezed thermal state defined in Eq. ([38|) . and 
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as well as 
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Further expanding the exponential term kt included in (C6) into sum series, and using the 
generating function of single-variable Hermite polynomials, 
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which leads to 
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Then using the recurrence relation of H n (x), 
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Eq.(C9) becomes 
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Thus we complete the derivation of WF Eq. flllj) by combing Eqs. (C4) and (C12). 
APPENDIX D: Derivation of (f45l) 
Substituting Eas. ([57f> - ([55| into flUD, we have 
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where T = (1 — e 2k *, (go,gi,g2, 93) and (x, w, G, Ai, A2) are defined in Eqs. (08]) and 
(f4Uj) . respectively. In a similar way to deriving Eq. fllll) . we can further put Eq.(Dl) into 
Eqs.gSD-glD. 

APPENDIX E: Derivation of (l57h 

Then employing Eqs. (|38p and (C6) as well as the integration formula (B7), we can treat 
the integration in a similar way to deriving Eq. ljlip . 
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Similarly to deriving Eq.(Bll), we have 
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and B 2 = B\- given in Eq.flgS}, which leads to Eq.fl57|). 
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